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We present results for mesonic propagators in temporal and spatial direction and for topological properties at 
T below and above the deconfining transition in quenched QCD. We use anisotropic lattices and Wilson fermions. 



1. Questions concerning QCD at T > 

With increasing temperature we expect the 
physical picture promoted by QCD to change ac- 
cording to a phase transition. Questions related 
with this transition are: whether chiral symmetry 
restoration and deconfinement represent the same 
phenomenon or succeed by very near but different 
transitions; how can one characterize the change 
in the vacuum structure (topology, monopoles, 
vortices, etc); which are the precise changes in 
the properties of hadrons; etc. The present paper 
reports the results from a project going on over a 
number of years and analyzing some of these ques- 
tions. Although this analysis is not closed, the 
present stage achieves a certain description which 
we shall present here together with the physical 
picture suggested by it and further problems. 

The primary question in this analysis has been 
the structure of the hadron correlators and its 
dependence on temperature. A second question 
raised in the later phase of the project concerns 
the relation to the topological properties. Both 
these problems have been treated in quenched 
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QCD. Both these questions will be pursued to 
completion. We also intend to extend preliminary 
studies for full QCD to a systematic analysis. 

Since T > breaks Lorentz symmetry the 
euclidean formulation of the finite temperature 
problem specifies a "temperature" (euclidean 
time) axis (see, e.g, Q). Therefore we need to 
investigate hadronic correlators with full "space- 
time" structure, in particular the propagation in 
the euclidean time direction. The latter, however, 
puts special problems because of the inherently 
limited physical length of the lattice, l T = ij,. 
To preserve at least a reasonably fine discretiza- 
tion of the time axisQ (needed in order to follow 
the details of the time-dependence of the corre- 
lators) we chose to use anisotropic lattices with 
a<r/«T = £ > 1 which can ensure the above re- 
quirement without increasing too much the to- 
tal volume of the lattice (section 2). To separate 
ground and excited states we need a more refined 
definition of the hadronic propagators. The strat- 
egy used will be described in section 3. Sections 
4 and 5 present our main results for the hadronic 
correlators and section 6 discusses the tempera- 
ture dependence of the topological properties of 
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the configurations. Section 7 is reserved for con- 
clusions and outlook. 

2. T > QCD on anisotropic lattices 

Anisotropic lattices are realized by introducing 
different space-space and space-time couplings, 
e.g. for the Wilson Yang-Mills action: 

Sym = -| (-ReTrOn + 1 Re r TiU a ^\ (1) 

This produces a cut-off anisotropy £ = 777(^,7), 
the temperature being given by T = j- = w ^ . 

The function ^(/3, 7) depends on the dynamics. 
To determine it correlation functions in different 
directions are calculated for given /?, 7 at T = 
and required to show isotropy under rescaling of 
the time distances by £ ("calibration"): 



I%(z) = FZ(t = tz) 



(2) 



Once the Yang-Mills calibration has been per- 
formed, we calculate hadron correlators from the 
action (X^ are lattice translation operators): 



S F = 2k (T #W#, W = 1 



V {^T+U l T l + lF T+U i T^j + "h.c", 
k^ 1 = 2(rn + 3 + 7 F ) = + j F - 1, 



(3) 
(4) 
(5) 



r± = i± 7|t , 7 " = i 



and tune 7^ to obtain the same £ when the phys- 
ical isotropy condition eq. (|^) is applied to these 
correlators. If the action leads to strong artifacts, 
(^) cannot be fulfilled simultaneously for all ob- 
servables and the result depends in particular on 
the action used. In [gj we compared the calibra- 
tion for various actions, see also jjj. 

Our quenched QCD analysis uses Wilson ac- 
tion with f3 = 5.68 and 7 = 4 on lattices of 
12 3 x N T with N T = 72, 20, 16 and 12. 
We have a a ~ 0.24 fm. The critical tempera- 
ture is found from the behavior of the Polyakov 
loop susceptibility to be at N T slightly above 18, 
which gives for the above lattices the tempera- 
tures T ~ 0, 0.93T C , 1.15T C and 1.5T C (see [|). 
The results to be reported below correspond to 
two sets of analysis: 
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Figure 1. Fermionic calibration. Three values of 
7f are tried with fixed K a = 0.081, and the value 
7f = 4.05 is adopted. 

- Set- A uses Wilson loops to perform the Yang- 
Mills calibration. We find a value for £ ranging 
between 5.3 for wide Wilson loops in eq. (||) and 
6.3 for thin Wilson loops. For the analysis an 
average £ = 5.9 is used and we find 7^ = 5.4. 
Because of the ambiguity in £ we use this set 
only to investigate the temperature dependence of 
the time-correlators but did not attempt to make 
a systematic comparison with space (or "screen- 
ing") correlators. 

- Set-B corresponds to a calibration performed for 
the static quark potential, which is a more ade- 
quate physical observable for requiring isotropy 
than the Wilson loops themselves ||. The am- 
biguity in fixing £ is reduced, we find a £ ~ 5.3 
and we can compare correlators corresponding to 
both time and space propagation. 

In Fig. [l] we present the calibration for Set-B. 
The various parameters are given in the table. 

3. Strategy 

Increasing the temperature is expected to in- 
duce significant changes in the structure of the 
hadrons. In particular, the mesonic peak in the 
spectral function of the correlators may move and 
develop a width. The question is, how to observe 
this change by distinguishing it from the influence 
of higher excitations which will not be dumped 
enough at the time scale of the inverse tempera- 
ture. 
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Two pictures are frequently used to describe 
the intermediate and the high temperature 
regimes: the weakly interacting meson gas and 
the quark-gluon plasma. In the first regime (as- 
suming for simplicity that we only deal with pi- 
ons) we obtain for the euclidean pion propagator 
at finite inverse temperature /? (t > 0): 

G^-coshOa^-t)), (6) 

with E n the energy of the n— pion state and 
A n = E n + E n — E n+ i measuring the interac- 
tion. As long as this interaction is weak we can 
simulate the changes by a shift and possibly a 
widening of the peak in the spectral function. If 
the changes are large - as they should be if the 
quark-gluon plasma regime becomes dominant - 
the above description breaks down. 

These genuine temperature effects should now 
be distinguished from the admixture of excited 
states in the mesonic channel under considera- 
tion. Our strategy is the following: we try to fix 
at zero temperature a mesonic source which gives 
a high projection onto the ground state. Then 
we use this source to determine the changes in- 
duced by the temperature on the ground state 
so defined. From the above considerations this is 
a reasonable procedure in the frame of a weakly 
interacting meson gas and it is justified as long 
as the observed changes are not too large. The 
appearance of large changes will then signal the 
breakdown of the weakly interacting gas picture 
and there we shall try to compare the observa- 
tions with those given by the quark-gluon plasma 
picture. We stress that in this case we do no 
longer have a good justification to use that source 
as representative of the pion. Our conjecture is 
that it still projects onto the dominant low en- 
ergy structure in the spectral function but further 



analysis are necessary to test this conjecture. In 
all the figures only statistical errors will be dis- 
played since we have no estimation for such sys- 
tematic effects. 

The correlators investigated are of the form: 

z yi,y2 

(Tr [7 5 5( yi , 0; z, t) l5 S(y 2 , 0; z + x, t)]) (8) 

For practical reasons we do not construct a 
genuine mesonic source but use smeared quark 
sources with point and exponential ansatze: 

w(yi,y 2 ) = wi(yi)w 2 (y 2 ) (9) 

w(y) ~ S(y) (point), w(y) ~ exp(-ay p ) (exp.) 

For the exponential source we fix the parameters 
a, p from the wave function (the dependence on x) 
measured with point source at T = 0. The results 
of a variational analysis using point-point, point- 
exponential and exponential-exponential quark 
sources indicate that the latter ansatz projects 
practically entirely on the ground state. This is 
well seen from the effective mass plots given in 
Fig. ^|. Therefore we use for our investigations 
throughout the exponential-exponential source. 

4. Temperature dependence of the time 
correlators 

We first present results for the effective mass of 
the Set- A data obtained from the it and p time- 
propagators at T ~ 0.93T C , 1.15T C and 1.5T C - see 
Fig. |[ Comparing with the results at T = (Fig. 
||) we notice practically no change at T ~ 0.93T C 
and relatively little change at 1.15T C . In con- 
tradistinction, the 1.5T C plot shows rather large 
effects: not only does the effective mass depend 
strongly on t and becomes significantly larger, but 
the 7r and p reverse their positions. 
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Figure 2. Effective pion mass vs t for various 
sources at T = 0: point-point, point-exponential 
and exponential-exponential. 
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Figure 4. The pion "wave functions" vs r for 
n a = 0.081 at N t = 72 and N t = 12. They are 
measured both in (0,0,1) and in (1,1,1) directions. 
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Figure 3. Effective it and p mass at T ~ 0.93, 
1.15 and 1.5T C (open symbols) vs t. Also shown 
are the effective masses from the same correlators 
calculated using free quarks. 

Of course, one should only compare distances 
up to 6, or about 1.4 fm, the largest distance avail- 
able at 1.5T C . Performing a fit to the propagators 
in this region and using a resonance ansatz shows 
a width developing above T c , leading to a very 
flat structure at 1.5T c f2] 

Since at 1.5T C there is every indication that the 
weakly interacting pion gas picture is lost we tried 



3 It is well known that due to the KMS condition poles 
cannot be moved off the real axis in an arbitrary way - 
see 111, see also |5|. It is possible, however, to write down 
an ansatz for a state with non- vanishing width which does 
not violate the KMS condition 



to see to what extent the situation is now com- 
patible with the quark-gluon picture. For this 
we compare the mesonic correlators with correla- 
tors obtained by using, in the same channels and 
with the same source, free quark propagators: Sq 
instead of S in (||) - see eq. (5) in (notice 
that this equation should be corrected at t = 
by adding a term (tanh(r?AT/2) + 7 4 )/2A). The 
results are also shown in Fig. [| The compar- 
ison cannot be made quantitative enough, since 
we cannot exactly determine what quark mass 
we should use for the free quarks. The data in 
the figure correspond to massless quarks, higher 
masses do not change the general aspect. The 
conjecture suggested by the comparison is that 
at 1.5T C a strong quark-gluon component has de- 
veloped and the mesons are in equilibrium with 
their decay products, quarks and gluons. 

5. Screening masses and quark mass de- 
pendence 

Further analysis of the pole and screening 
masses, and their dependence on the quark mass 
are carried out on the Set-B, on which three val- 
ues of quark mass and higher statistics (30 con- 
figurations) are available. 

The pole masses m = are extracted from 
the correlators in time (temperature) direction 
using again smeared quark propagators with the 
exp. - exp. source tuned at T — 0, as described 
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in section 3. We stress again that at high tem- 
perature there is a certain arbitrariness in con- 
sidering the state on which this source projects 
as the dominant low energy state. The effective 
masses show a similar behavior with those from 
Set- A. For further considerations we fit the cor- 
relators to single hyperbolic functions in the re- 
gions: 27-45 (T = 0), 7-13 (T ~ 0.93T C ), 6-10 
(T ~ 1.15T C ) and 4-8 (T ~ 1.5T C ), which are cho- 
sen observing the effective mass plots. The results 
of K a = 0.081, our heaviest quark mass, are shown 
in Fig. |]and Fig. |5|. In Fig. |we show the depen- 
dence of the correlators on the distance between 
quark and anti-quark at the sink at T = and at 
1.5T C . The wave function is seen to become more 
flat at higher temperature. 

The screening masses are extracted from 
the correlators in z-direction. At T — 0, = 
to ( - T ) is ensured by tuning 7^. At finite tempera- 
ture, in general the screening mass and the pole 
mass do not coincide. Since the physical distance 
in the space direction is reasonably large we ex- 
tract the screening masses from correlators with 
point-point source in the range 5-7 at all temper- 
atures (although at T — no clear plateau is yet 
observed, the result agrees within a few percents 
with the mass obtained with a Wuppertal source). 
At T ~ 0.93T C the z-correlators are very similar 
to those at T = 0. Above T c , the effective masses 
seem to reach the plateau beyond z — 4. The 
results for m (<7) at n a = 0.081 are also shown in 
Fig. |. 

Both pole and screening masses are extrapo- 
lated to the chiral limit from the 3 quark masses 
analyzed. The extrapolation is done using k of 
eq. (Q). Below the critical temperature the pion 
masses squared are extrapolated, in all other 
cases, the masses themselves are extrapolated lin- 
early in 1/k. The extrapolation of the pole masses 
at T = gives the critical value n c — 0.17138(20) 
and a^ 1 = 1.18 GeV from the p mass using 
£ = 5.3. The result of the screening mass extrap- 
olation is almost the same - see Fig |[ The large 
difference between a~ 1 {p) and a~ l from the heavy 
quark potential (see the table) indicates signif- 
icant discretization effects for our lattices as ex- 
pected. The results are summarized in the second 
plot in Fig. ||. We find large differences between 
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Figure 5. Temperature dependence of the pole 
and the screening masses for n a = 0.081(top fig- 
ure), and at the chiral limit (bottom). The verti- 
cal gray lines roughly represent the critical tem- 
perature. 

the pole and the screening masses above T c both 
at finite quark masses and in the chiral limit. We 
should like to point out that a similar behavior 
is obtained in an effective model approach with 
NJL Lagrangian ||. 

6. Topological properties 

To determine the topological properties of the 
configurations we use improved cooling - see 
||. Unlike usual cooling with Wilson action this 
method preserves instantons at scales larger than 
a certain dislocation threshold po- F° r anisotropic 
lattices the improved action is modified accord- 
ingly ensuring isotropic distributions for the sizes 
and distances of instantons at T = 0. From ear- 
lier investigations we expect a scale of about 0.4 
- 0.7 fm to be relevant for the topology of 5^7(3). 
With a dislocation threshold po ~ 2a ~ 0.5 fm 
the lattices used here are too coarse for a de- 
tailed analysis. Nevertheless we obtain reason- 
able data for the susceptibility in agreement with 
Witten-Veneziano formula below T c and dropping 
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Figure 6. Chiral extrapolation of the pion 
screening masses at various T. 

to zero at high temperature, in agreement with 
other analysis on finer lattices. An interesting 
quantity is the ratio R(k) = Q^{kV)/Q 2 s (kV) be- 
tween the charge in subvolumes kV obtained by 
slicing the lattice in either a spatial or a tempo- 
ral direction. The temperature effect shows up in 
an anisotropy between temporal and spatial slic- 
ing which sets in when the dominant sizes begin 
to feel the limited temporal extension. Just be- 
low T c R(l/2) in the charge zero sector is about 
0.7, while above T c (where all configurations have 
Q = 0) #(1/2) grows fast to about 3. This be- 
havior is compatible with a random distribution 
of instantons of rather large sizes below T c and 
a tendency to form pairs above T c . Generally, 
there seems to be a qualitative difference between 
the instanton population below and above T c : at 
similar densities of pairs the instantons of the low 
temperature are reasonably well fitted by the 't 
Hooft ansatz, while at high temperature we seem 
to deal with topological charge fluctuations which 
hardly can accommodate the continuum picture 
of instantons. A more detailed account of these 
features will be given elsewhere. 

7. Conclusions and outlook 

The present analysis has produced results 
about the temperature effects on the vacuum 
and hadronic properties in quenched QCD which 
are consistent with the following picture: be- 
low T c the changes are small and gradual while 
above T c the changes increase strongly (but not 
abruptly) with the temperature; at temperatures 



of about 1.5T C the mesons have become unsta- 
ble and are strongly interacting with a signifi- 
cant quark-gluon plasma component; the changes 
in the vacuum structure follow a similar pattern 
and there is indication of close range opposite 
charge correlation (although there are no dynam- 
ical fermions). 

This picture needs to be further tested in or- 
der to remove the uncertainties still affecting 
the present analysis. This concerns particularly 
the question of an unambiguous definition of the 
hadron states at high temperature. In the fur- 
ther developments we shall also try to extract di- 
rectly information about the spectral function of 
the hadronic correlators fLOfl , use larger lattices 
such that we can both have larger physical dis- 
tances and a smaller lattice spacing (in spatial 
directions) for more refined studies of the topol- 
ogy, and aim, of course, to an improved statistics. 

The calculations have been done on API 000 
at Fujitsu Parallel Computing Research Facilities 
and Paragon at INSAM, Hiroshima Univ. 
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